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This paper presents an under-appreciated way to conceptualize stationary black holes, which we 
call the river model. The river model is mathematically sound, yet simple enough that the basic 
picture can be understood by non-experts. In the river model, space itself flows like a river through 
a flat background, while objects move through the river according to the rules of special relativity. 
In a spherical black hole, the river of space falls into the black hole at the Newtonian escape velocity, 
hitting the speed of light at the horizon. Inside the horizon, the river flows inward faster than light, 
carrying everything with it. We show that the river model works also for rotating (Kerr-Newman) 
black holes, though with a surprising twist. As in the spherical case, the river of space can be 
regarded as moving through a flat background. However, the river does not spiral inward, as one 
might have anticipated, but rather falls inward with no azimuthal swirl at all. Instead, the river 
has at each point not only a velocity but also a rotation, or twist. That is, the river has a Lorentz 
structure, characterized by six numbers (velocity and rotation), not just three (velocity). As an 
object moves through the river, it changes its velocity and rotation in response to tidal changes in 
the velocity and twist of the river along its path. An explicit expression is given for the river field, 
a six-component bivector field that encodes the velocity and twist of the river at each point, and 
that encapsulates all the properties of a stationary rotating black hole. 

PACS numbers: 04.20.-q 



I. INTRODUCTION 

As first pointed out in 1921 by Allvar Gullstrand 1 
and Paul Painleve^, the Schwarzschild^ metric can be 
expressed in the form 

ds 2 = -dt 2 s + (dr + f3dtff) 2 + r 2 (d9 2 + sm 2 6d4> 2 ) (1) 

where (3 is the Newtonian escape velocity, in units of the 
speed of light, at radius r from a spherical object of mass 
M 



and iff is the proper time experienced by an object that 
free falls radially inward from zero velocity at infinity. 

Although Gullstrand's paper was published in 1922, 
after Painleve's, it appears that Gullstrand's work has 
priority. Gullstrand's paper was dated 25 May 1921, 
whereas Painleve's is a write up of a presentation to 
the Academie des Sciences in Paris on 24 October 1921. 
Moreover, Gullstrand seems to have had a better grasp 
of what he had discovered than Painlcvc, for Gullstrand 
recognized that observables such as the redshift of light 
from the Sun are unaffected by the choice of coordinates 
in the Schwarzschild geometry, whereas Painleve, noting 
that the spatial metric was flat at constant free-fall time, 
dt ff — 0, concluded in his final sentence that, as regards 
the redshift of light and such, "e'est pure imagination de 
pretendre tirer du ds 2 des consequences de cette nature" . 

As shown in ||nj the Gullstrand-Painleve metric 
provides a delightfully simple conceptual picture of the 
Schwarzschild geometry: it looks like ordinary flat space, 
with the distinctive feature that space itself is flowing 




FIG. 1: (Color online) The fish upstream can make way 
against the current, but the fish downstream is swept to 
the bottom of the waterfall. Figure 1 oJ- presents a similar 
depiction. 



radially inwards at the Newtonian escape velocity. The 
place where the infall velocity hits the speed of light, 
/3 = 1, marks the horizon, the Schwarzschild radius. 
Inside the horizon, the infall velocity exceeds the speed 
of light, carrying everything with it. 

Picture space as flowing like a river into the 
Schwarzschild black hole. Imagine light rays, photons, 
as fishes swimming fiercely in the current. Outside the 
horizon, photon-fishes swimming upstream can make way 
against the flow. But inside the horizon, the space river 
is flowing inward so fast that it beats all fishes, carrying 
them inevitably towards their ultimate fate, the central 
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singularity. 

The river model of black holes offers a mental picture 
of black holes that can be understood by non-experts 
(at least in the spherical case) without the benefit of 
mathematics. It explains why light cannot escape from 
inside the horizon, and why no star can come to rest 
within the horizon. It explains how an extended object 
will be stretched radially by the inward acceleration of 
the river, and compressed transversely by the spherical 
convergence of the flow. It explains why an object 
that falls through the horizon appears to an outsider 
redshifted and frozen at the horizon: as the object 
approaches the horizon, light emitted by it takes an ever 
increasing time to forge against the inrushing torrent of 
space and eventually to reach an outside observer. The 
river model paints a picture that is radically different 
from the Newtonian picture envisaged by Michell (1784) 6 
and Laplace (1799) 7 .. 

The picture of space falling like a river into a black 
hole may seem discomfortingly concrete, but the aetherial 
overtones are no more substantial than in the familiar 
cosmological picture of space expanding (see e.g. p. 237 
of Greene 20042). 

As reviewed by Visser (1998, 2003) 9 ' 10 and 
by Martel & Poisson (2001) 11 , the Gullstrand- 
Painleve metric has been discovered and rediscovered 
repeatedly 12 i 13 i 14 i 15 i 16 i 17 i 18 i 19 i 2Q i 21 i 22 i 23 . Surprisingly, 
the Gullstrand-Painleve metric is widely neglected in 
texts on General Relativity. An admirable exception is 
the text "Exploring Black Holes" by Taylor & Wheeler 
(20 00) 24 , which devotes an entire section, Project B, 
to the Gullstrand-Painleve metric, calling it the "rain 
frame" (the metric itself appears on page B-13). Taylor 
& Wheeler attribute (page B-26) the idea for the rain 
frame to the book by Thorne, Price & MacDonald 
(1986) 2 », page 22 and elsewhere, although the metric 
does not appear explicitly in the latter book. 

It has been recognized for decades that some aspects 
of general relativity can be conceptualized in terms 
of flows. In the ADM (1962) formalism^ 6 , (see 
e.giSl for a pedagogical review), one considers fiducial 
observers - FIDOs^ 5 - - whose worldlines are orthogonal 
to hypersurfaces of constant time. The shift vector in 
the ADM formalism is just the velocity of these FIDOs 
through the spatial coordinates. Alcubierre (1994) 28 
constructed his famous warp-drive metric by positing a 
superluminal (faster-than-light) shift vector. 

In a seminal, albeit initially unremarked paper, Unruh 
(1981)2£ (see 3 ^ for a comprehensive review) pointed out 
that the equations governing sound waves propagating in 
an inviscid, barotropic (pressure is a specified function 
of density), irrotational fluid are the same as those for 
a massless scalar field propagating in a certain general 
relativistic metric. Unruh showed that this implied that 
sound horizons would emit Hawking radiation in much 
the same way as event horizons in black holes, and he 
proposed that Hawking radiation might be detected from 
sonic black holes, or "dumb holes" , in the laboratory. 



As excellently reviewed by Barcelo et al. (2005)^, 
Unruh's paper in due course led to a now thriving 
industry on "analog gravity" , in which fluid flows with 
prescribed velocity fields simulate general relativistic 
spacetimes. The primary aim of the work on analog 
gravity is to try to understand, and perhaps in the 
not-too-distant future to probe experimentally, quantum 
gravity through sonic analogs. 

It is generally assumed that the fluid, or river, 
analogy applies to a limited class of general relativistic 
spacetimes, those in which the metric can be expressed 
up to an overall factor (a "conformal" factor) in terms of 
a shift vector (the velocity of the river) on an otherwise 
flat background space. The 3-dimensional shift vector 
and the conformal factor provide 4 degrees of freedom, 
whereas at least 6 degrees of freedom are required to 
specify an arbitrary spacetime (the metric has 10 degrees 
of freedom, of which 4 are removed by an arbitrary 
coordinate transformation). As a corollary, it has been 
thought that any general relativistic geometry admitting 
a fluid analog must necessarily be (up to a conformal 
factor) spatially flat at constant tim o 31 i 32 , as indeed is 
the case in the Gullstrand-Painleve metric. 

In particular, it has been thought that no river model 
for stationary rotating black holes exists^, since the 
Kerr-Newman geometry does not admit conformally flat 
slices 3 ^ 3 !. 

In the present paper we started from a somewhat 
different conceptual picture. We noticed that fishes 
swimming in the Gullstrand-Painleve river moved ac- 
cording to the rules of special relativity, being boosted by 
tidal differences in the river velocity from place to place. 
We wondered, might there be an analogous behavior for 
rotating black holes? It came as a magical surprise, flllll 
that the answer is yes, from this perspective there is a 
river model of the Kerr-Newman geometry. The rotating 
analog of the Gullstrand-Painleve metric proves to be 
(as expected 3 ^) the Doran (20 00) 35 form of the Kerr- 
Newman metric. The new feature that emerges from the 
mathematics is that the river of a rotating black hole is a 
fully 6-dimensional Lorentz river, with a twist as well as a 
velocity. Just as a velocity is a generator of a space-time 
rotation (a Lorentz boost), so also a twist is a generator 
of a space-space rotation (an ordinary spatial rotation) . 
As a fish swims through the Doran river, it is not only 
boosted but also rotated by tidal differences in the river 
velocity and twist from place to place. 

This novel point of view leads to a different notion 
of what is meant by the flat background space through 
which the river flows and twists. Mathematically, 
the essential feature of the river model appears to 
be equation l|71|) . which states that the connection 
coefficients, expressed in locally inertial frames comoving 
with the infalling river of space, should equal the ordinary 
(non-covariant) gradient of the river field. 

The property that the tetrad connection coefficients 
are equal to the ordinary gradient, in Doran-Cartesian 
coordinates, of the river field, essentially defines what we 
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mean by the background space in the river model being 
flat. This feature appears to be a special property of 
stationary black holes. How this idea emerges from the 
mathematics is examined in Wl Fl and revisited in fllll II 
We emphasize that the background being flat does not 
mean that the metric is spatially flat, although the latter 
is also true in the case of spherical black holes. The 
notion that there is a sense in which stationary rotating 
black holes admit a flat background coordinatization 
might have application to numerical general relativity, 
for example in setting up initial conditions containing 
rotating black holes, where traditional conformal imaging 
and puncture methods that assume a conformally flat 3- 
geometry are too restrictive to admit Kerr black holea 3 ^. 

Throughout this paper we adopt the sign conventions 
and ordering of indices of Misner, Thorne and Wheeler 37 . 



II. SPHERICAL BLACK HOLES 

In this section we consider spherically symmetric black 
holes, and we justify the assertion that the Gullstrand- 
Painleve metric, equation JQl, can be interpreted as 
representing a river of space falling radially inward at 
velocity (3. We demonstrate two features that are the 
essence of the river model for spherical black holes: 
first, that the river of space can be regarded as moving 
in Galilean fashion through a flat Galilean background 
space [eqs. (THTl) and <fT5|) ]. and second, that as a freely- 
falling object moves through the flowing river of space, 
its 4- velocity, or more generally any 4-vector attached to 
the freely-falling object, can be regarded as evolving by 
a series of infinitesimal Lorentz boosts induced by the 
change in the velocity of the river from one place to the 
next [eq. (|f 8f) ]. Because the river moves in a Galilean 
fashion, it can, and inside the horizon does, move faster 
than light through the background. However, objects 
moving in the river move according to the rules of special 
relativity, and so cannot move faster than light through 
the river. 



A. Mathematics of the river model 

In general, a spherically symmetric metric of the form 
(units c = G = 1) 



ds 2 



[l-2M(r)/r]dt 2 + 



dr 2 



[1 - 2M(r)/r] 
+ r 2 (d9 2 +sm 2 0d(f> 2 ) (3) 



can be expressed in the Gullstrand-Painleve form Q 
with infall velocity 



/3(r) 



2M(r) 



1/2 



(4) 



the free-fall time tg being 



t a = t 



1-I3 2 



■ dr 



(5) 



The velocity (3 is commonly called the shift in the ADM 
formalismSiiS 7 ., but in this paper we refer to (3 as the 
river velocity. The river velocity (3 is positive for a black 
hole (infalling), negative for a white hole (outfalling). 
Horizons occur where the river velocity (3 equals the 
speed of light, 



P = ±1 



(6) 



with (3 = 1 for black hole horizons, and (3 = — I for 
white hole horizons. The Reissner- Nordstrom metric for 
a spherically symmetric black hole of mass M and charge 
Q takes the form @ with mass M(r) interior to r, the 
so-called Misner-Sharp mass (Misner & Sharp 1964 38 ; 
Weinberg 1972, p 30C 3 £), given by 



O 2 

M(r) = M -— . 
v ' 2r 



(7) 



However, the river velocity [3 can also be considered to 
be a more general function of radius r. In All Bl we will 
return briefly to the Reissner-Nordstrom solution to see 
what its river looks like. 

To make the argument plainer, rewrite the Gullstrand- 
Painleve metric 
(a; , x 1 , x 2 
nates: 



in Cartesian coordinates x p 
! ) = {ts,x,y,z) instead of spherical coordi- 



ds 2 = v^idx^ ~ f3 fi dt s )(dx' J - I3 v dt s ) 
where rj^ v is the Minkowski metric, and 

(3" = (3 ( 



x 
r 



y 

r 



(8) 



(9) 



are the components of the radial river velocity. 

Let denote the basis of tangent vectors in the 
Gullstrand-Painleve-Cartesian coordinate system x^. By 
definition, the scalar products of the tangent vectors 
constitute the metric 



(10) 



Let = dx^/dr denote the 4- velocity of a particle 
falling freely (not necessarily radially) in the geometry, 
r being the proper time experienced by the particle. 
In particular, observers who free-fall radially from zero 
velocity at infinity have 4-velocity 



v£ = (l,f3\[3 2 ,{3 3 ) 



fill 



Such observers are comoving with the inflowing river 
of space. Let j m , and associated local coordinates 
gm _ (^°,^ 1 ,^ 2 ,^ 3 ), denote a system of locally inertial 
orthonormal frames, tetrads, attached to observers who 
free-fall radially from zero velocity at infinity. Here 
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and throughout this paper we use latin indices to 
signify tetrad frames, reserving greek indices for curved 
spacetime frames. Orthonormal means that the scalar 
products of the tetrad basis at each point of spacetime 
form the Minkowski metric 



~fm ' in Vn 



(12) 



That the tetrad frames move with the radially free-falling 
observers without precessing requires that the vectors 
7m be 'parallel-transported' along the worldlines of these 
observers, that is, 



a d"y m 
ff dx^ 



(13) 



Assume without loss of generality that the tetrad frames 
7 m are aligned with the Gullstrand-Painleve-Cartesian 
frame at infinity. Then the tetrad frames 7 m are related 
to the Gullstrand-Painleve basis at each point by 



7o 

7* : 



: 9o 
9i 



(i = 1,2,3) 



(14) 



which is most easily deduced from equations (|10fl and 
(|12[) and considerations of symmetry, the result being 
confirmed by checking that equation 1)13(1 is then true. 
Remarkably, the relations (|14fl are those of a Galilean 
transformation, which shifts the time axis by velocity /3 
along the direction of motion, but leaves unchanged both 
the time component of the time axis and all of the spatial 
axes. 

The 4-velocity u m of a freely-falling particle with 
respect to the the tetrad frame 7 m at the position of 
the particle follows from u rn ~f m — v^g^, which implies 



u l = v l -P l v° (i = 1,2,3) 



(15) 



Physically, the 4-velocity u rn is the 4-velocity of the 
particle relative to the inflowing river of space. For 
example, the spatial components u l of the 4-velocity are 
zero if the particle is becalmed in the river. Again, the re- 
lations (|15|) resemble those of a Galilean transformation, 
which shift only the spatial components of the vector, 
while leaving the time component unchanged. The only 
non-Galilean (relativistic) feature of equations (f 1 51) is 
that the 4-velocities u rn and u' 1 are derivatives with 
respect to proper time. But proper time is a property 
of the objects moving in the river, not of the river itself. 
Objects moving in the river move through it according 
to the rules of special relativity. The river itself flows in 
Galilean fashion through a flat Galilean background. 

We have demonstrated, equations (|14|) and H15|l . the 
first of the two claimed features of the river model for 
spherical black holes, that the river of space moves in 
Galilean fashion through a flat Galilean background. We 
proceed to demonstrate the second feature of the river 
model, that objects moving in the river of space move 
according to the rules of special relativity, being Lorentz 



boosted by tidal differences in the river velocity from 
place to place. 

The tetrad frames have been constructed so that an 
observer who free-falls from zero velocity at infinity 
finds their own frame aligned at all times with the 
tetrad frame. But in general another observer who 
free-falls along a different geodesic will find their 
own locally inertial frame becoming misaligned with 
the tetrad frame. This misalignment is determined 
mathematically by the equations of motion of objects, 4- 
vectors, expressed with respect to the tetrad frame. Let 
p = p m -fm = p^gfj, be a 4-vector. Its components p m in 
the tetrad frame are related to those p^ 1 in the coordinate 
frame by p m = 5™p^ — (3 m p°. For a 4-vector in free-fall, 
the equations of motion for the components p m in the 
tetrad frame are (see ijlll Efl 



dp k 



T« mn u n p m = 



(16) 



dp 




dT 




dp 1 


dP l j o 


~d7 


~dxT Up 



where T mrl are the tetrad frame connection coefficients, 
the tetrad frame analog of the coordinate frame Christof- 
fel symbols. In the present case of spherically symmetric 
black holes, the non-zero tetrad frame connection 
coefficients prove to be given by the spatial gradient of 
the river velocity (see Sjlll Hp 

8 ft* 

1" (i,3 = 1,2,3). (17) 

From equations (|16fl and (|17|) it follows that 



(18) 



The summations over paired indices in equations l|18|) 
are formally over all four indices 0, 1, 2, 3, but in practice 
reduce to sums over just the three spatial indices 1,2,3 
since, first, the infall velocity has zero time component, 
ft = 0, as we have defined it, and second, the infall 
velocity has zero time derivative, dB^/dx = 0. 

In the context of the river model, the equations of 
motion (|18fl have the following interpretation. In an 
interval St of proper time, a particle moves a distance 
Sx l — v 1 5t in the background Gullstrand-Painleve- 
Cartesian coordinates, and a proper distance <5£* = 
u 1 5t = Sx l — (3 l Stff relative to the infalling river of space. 
The proper distance <5£ l equals the distance 5x l minus 
the distance f3 l Ste moved by the river. In Gullstrand- 
Painleve-Cartesian coordinates, the velocity B l of the 
infalling river at the new position differs from the velocity 
at the old position by 8x 3 d(3 l /dx J However, in the 
river model, a particle moving in the river sees not the 
full change in river velocity relative to the background 
coordinates, but only the tidal change 



5P" 



8B l 
dxi 



6& 



(19) 
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in the river velocity relative to the infalling locally inertial 
river frame. For example, if the particle is comoving with 
the inflowing river, so that <5£* = 0, then the particle sees 
no change at all in the river velocity as time goes by, 
5/3' 1 = 0. The infinitesimal tidal change 5f3 l in the river 
velocity induces a Lorentz boost in the 4- vector p m 



P 
P* 



P°-6f3 lP * 
Sf3 l P ° 



P 



(t = l,2,3) 



(20) 



Equations (|19|l and (|20|l reproduce the equations of 
motion H18|l. 

We have thus demonstrated the second of the claimed 
features of the river model for spherical black holes, that 
as a particle moves through the river of space, its 4- 
velocity, or more generally any 4- vector attached to it, 
is Lorentz boosted by tidal changes in the river velocity 
along its path. 



B. Reissner-Nordstrom metric 

We conclude this section by commenting on the 
Reissner-Nordstrom (RN) metric for a spherical black 
hole of mass M and charge Q. In this case the mass 
M(r) = M — Q 2 /(2r) interior to r, equation (7J, can 
be interpreted as the mass M at infinity, less the mass 
(E 2 /8tt) 4nr 2 dr = Q 2 /(2r) contained in the electric 
field E = Q/r 2 outside r. The RN geometry exhibits 
both outer and inner horizons r + and r_ 



r ± =M±(M 2 -Q 



2U/2 



(21) 



The inflow velocity (3 hits the speed of light at the outer 
horizon r+, reaches a maximum velocity between the 
outer and inner horizons, slows back down to the speed 
of light at the inner horizon r_ , and slows all the way to 
zero velocity at the turnaround radius 



To 



2M 



(22) 



At this point the flow of space turns around, accelerates 
back outward through another inner horizon, the so- 
called Cauchy horizon, into a white hole, and bursts 
through the outer horizon of the white hole into a new 
universe. 

Sadly, the RN solution is not realistic, and its promise 
of passages to other universes is moot. The RN solution 
describes the geometry of empty space surrounding 
a classical point charge of infinitesimal size, and as 
such consists of a gravitationally repulsive singularity 
of infinite negative mass sheathed in an electric field 
containing an infinite positive mass balanced so as to 
yield a finite mass M at infinity. See^S for an entry to 
the literature on the intriguing subject of what really 
happens inside charged black holes. 

The infall velocity /3 is imaginary inside the turnaround 
radius rp of the Reissner-Nordstrom geometry, the 



interior mass M(r) being negative inside this radius. 
This might be considered a defect of the river model, 
but it might also be considered an asset, signalling the 
presence of unphysical negative mass. Whatever the case, 
the mathematical formalism remains valid even where the 
river velocity (3 is imaginary. 



III. ROTATING BLACK HOLES 

Does the river model work also for stationary rotating 
black holes? As will be shown in this section, the answer 
is yes: there is a river of space, and it moves through a 
flat background, and fishes move through the river special 
relativistically as though they were being carried with 
it. But the river has a surprising twist. One might have 
anticipated that the river would spiral into the black hole 
like a whirlpool, but that is not the case. Rather, the 
river velocity has no azimuthal component at all. Instead 
of a spiral, the river possesses, besides a velocity at each 
point, a rotation, or twist, at each point. The river is 
characterized not by three numbers, a velocity vector, 
but by six numbers, a velocity vector and a twist vector. 
As a fish swims through the river, it is Lorentz boosted by 
gradients in the velocity of the river, and rotated spatially 
by gradients in the twist of the river. 

A key result of this section is the expression (|72|l for 
the river field oJkm- This is a 6-component bivecto»2Ziii 
field, antisymmetric in its indices km, whose electric 
part specifies the river velocity, and whose magnetic part 
specifies the river twist. The river field Ukm encapsulates 
all the properties of a stationary rotating black hole. 

How can a river move and twist without spiralling? 
The answer to this conundrum is that, unlike the 
Gullstrand-Painleve case, the spatial metric is not flat, 
but sheared (Fig.|21). One can regard the twist in the river 
as inducing the shear in the spatial metric; or equally well 
one can regard the shear in the spatial metric as requiring 
a twist in the river. Whatever the case, the twist and the 
shear act together in just such a way as to ensure that 
locally inertial frames moving through the infalling river 
comove with the geodesic motion of points at rest in a 
small neighborhood of the frame (Fig.^J. 

Recall from special relativity that Lorentz transfor- 
mations are generated by a combination of changes 
in velocity, or Lorentz boosts, and spatial rotations. 
Lorentz boosts are rotations in a plane defined by a space 
axis and a time axis, while spatial rotations are rotations 
in a plane defined by two spatial axes. Gradients in the 
velocity of the river make the metric non-flat with respect 
to the time components, while leaving the spatial metric 
at constant time flat. Gradients in the rotation, or twist, 
of the river make the metric non-flat with respect to the 
spatial components, while leaving the time part of the 
metric flat, that is, the metric becomes — dt 2 + gijdx 1 dx J 
where g^ is a purely spatial metric. We see that the 
reason that the Gullstrand-Painleve metric for spherical 
black holes is flat along hypersurfaces of constant free-fall 
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time is attributable to the fact that the river has no twist 
component. However, the Gullstrand-Painleve river does 
have a velocity component, so the Gullstrand-Painleve 
metric is not flat in the time direction. For rotating black 
holes, the river has both velocity and twist components, 
and the metric is flat neither in time nor in space. 



A. Doran metric 

Proceed to the mathematics. Doran (2000)^ has 
pointed out that the Kerr-Newman metric for a rotating 
black hole of angular momentum a per unit mass (for 
positive a, the black hole rotates right-handedly about 
its axis) can be cast in oblate spheroidal coordinates 
(iff, r, 9, (be) in the form 



ds 2 



■dtl 



-\- p 2 d9 2 ~T n Mil 1/ m 

where (3{r) is the river velocity, 

R ee (r 2 + a 2 ) 1 ' 2 , p ee (r 2 + a 2 cos 2 I 



pdr (3R . o 

— — -| (dte — asm 9 dq>e , 

R p 

- R 2 sm 2 9 r^ 2 



a/2 



(23) 



(24) 



and the free-fall time ts and free-fall azimuthal angle 
(be are related to the usual Boyer-Lindquist time t and 
azimuthal angle (b by 



t ff = t 



be = 4> - a 



(idr 
l-(3 2 



Pdr 



R 2 (l-[3 2 ) 



(25) 



(26) 



As before, we adopt the convention that the river velocity 
(3 is positive for a black hole (infalling), negative for a 
white hole (outfalling) . Horizons occur (see AIII Dll where 
the river velocity (3 equals the speed of light 



= ±1 



(27) 



with (3 = 1 for black hole horizons, and (3 = — 1 for 
white hole horizons. Ergospheres occur where ds 2 = at 
dr = d9 = dtps = 0, which happens at 



' R 



(28) 



again with (3 — p/R for black hole ergospheres, and (3 — 
—p/R for white hole ergospheres. For a Kerr-Newman 
black hole with mass M and charge Q, the river velocity 
(3 is 



0(r) - 



(2Mr - Q 2 ) 1 / 2 



R 



(29) 



but for the present purpose the river velocity can be 
considered to be a more general function of the radial 



coordinate r. Note that the river velocity (3 defined here 
differs from Doran'sS* velocity by a factor of p/R- Doran 
defines the velocity to equal the magnitude (/^/J^) 1 / 2 = 
(3R/ p of the velocity vector (3^ given by equation lETT|) 
below, a seemingly natural choice. The point of the 
convention adopted here is that (3(r) is any and only a 
function of r, rather than depending also on 9 through 
p ee (r 2 + a 2 cos 2 9) 1 / 2 . Moreover, with the convention 
here, the river velocity is plus or minus one at horizons, 
equation l(77|) . as will be demonstrated below, mil Dl 

If the river velocity (3 is zero, then the metric (|23f) 
reduces to the flat space metric in oblate spheroidal 
coordinates. However, unlike the spherical case, the 
metric is not flat along hypersurfaces of constant free- 
fall time, dig = 0. 



B. Doran-Cartesian metric 

The Doran coordinate system turns out, S flHI Fl and 
IIII II to provide the coordinates of the flat background 
through which the river of space flows into the black 
hole. We therefore express the Doran metric in Cartesian 
coordinates x M = (x°, x 1 , x 2 , x 3 ) = (tg,x,y,z) — 
(iff , R sin 9 cos (be ,Rsm9 sin (be , r cos 9) with the rotation 
axis taken along the z-dircction: 

ds 2 = rj^idx" - f3 fl a K dx K )(dx l/ - f3 u a x dx x ) . (30) 

Here the components of the river velocity (3^ are 

^ = f Q _W_ 

p \ ' Rp ' Rp ' rp 

and a^dx 11 — dte — asin 2 9 dtfie has components 



(31) 



a n = ( 



' ' R 2 ' R 2 ' 



(32) 



The vector a M is related to the 4-velocity of the horizon, 
equation 1)460. and we refer to it as the azimuthal vector, 
since its spatial components point in the (negative) 
azimuthal direction, in the direction opposite to the 
rotation of the black hole. The spheroidal radial 
coordinate r is given implicitly in terms of x, y, z by 



2„2 



2/2, 2, 2 2\ 2 

r (x + y + z — a ) — a z 







(33) 



River tetrad 



In modelling black holes as an inflowing river of 
space, it is natural to work in the orthonormal tetrad 
formalism. Let denote the basis of tangent vectors 
in the Doran-Cartesian coordinate system x^, and let 
7 m , and associated local coordinates £ m = (£° , i^ 1 , £ 2 , £ 3 ) , 
denote a system of locally inertial frames, tetrads, 
attached to observers who free-fall from zero velocity 
(with zero angular momentum) at infinity. Such freely- 
falling observers are comoving with the infalling river of 



7 



space. They fall along trajectories of constant 9 and c/>s, 
and have 4-velocities Vg = (1, /? 2 , /3 3 ) in the Doran- 
Cartesian coordinate system. The scalar products of the 
tangent vectors at each point constitute the metric 
g M „, equation IjlOl) . while the scalar products of the tetrad 
vectors 7 m at each point form the Minkowski metric, 
equation (|12fl . If the tetrad frames 7 m are assumed, 
without loss of generality, to be aligned with the tangent 
vectors g M at infinity, then the relation between 7 m and 
is, as previously noted by Doran (2000)^, 



7tf t 



9t a + fig* 
g x + a x /3 l gi 
9y + a y P l g* 
9z 



(34) 



which may be confirmed by checking that the scalar 
products of the 7 TO so constructed form the Minkowski 
metric, and that their derivatives vanish along the 
worldlines of observers who free-fall from zero velocity 
at infinity, v^d^f m / dx^ — 0. If horizontal radial and 
azimuthal axes are defined by (7_> , 7j ) = (cos tfis *) x + 
sin (t>sl?y,- sin <p s ~f x + cos </> ff j y ) and likewise (g^ , g T ) = 
(cos ff g x + sin <p s g y , - sin <p s g x + cos (f>sg y ), then 



7^ = 9^ 

a sin 9 
7T = 9T R~P 9i 



(35) 



Equations (|34|l and (|35|l show that the time axis 7t ff is 
shifted by velocity f3 l gi, similar to the spherical case, 
equation l|14|l . but in addition the azimuthal axis 71 
is shifted by —{asm9/R)(3 l gi. Figure illustrates the 
horizontal radial and azimuthal axes 7^ and 71 at several 
points in the equatorial plane of a Kerr black hole. The 
azimuthal axes 7^ are tilted radially, in accordance with 
equation l|35|l . reflecting the fact that the spatial metric 
is sheared. 

Equations l|34|l may be abbreviated j m = e m ^g^ where 
e m M is the vierbein 

e TO " = <% + a m ^ (36) 
with 5^ n a Kronecker delta. The inverse vierbein e m u is 



1 m 



(37) 



That the product of the vierbein and its inverse given 
by equations (|36(l and (|37|l is indeed the unit matrix, 
e m M e" M = <S™ and e" 1 (I e m " = (5^, follows from the 
orthogonality of the azimuthal and velocity vectors a M 
and /3 M , namely = 0. The vectors a m with a latin 

index in the vierbein l|3t)|l and (3 m with a latin index in 
the inverse vierbein l|37|l are defined by 



Sta b 



and transform with the tetrad frame -y m rather than the 
coordinate frame g^. The coordinates of a m and (3 m 
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FIG. 2: (Color online) Sets of horizontal radial and azimuthal 
tetrad axes 7^ and 71, equations 1351 . in the equatorial s- 
y plane ((9 = 7r/2) of an uncharged (Kerr) black hole with 
angular momentum per unit mass a = 0.96, plotted in Doran- 
Cartesian coordinates. The azimuthal axis at each point is 
tilted radially, reflecting the fact that the spatial metric is 
sheared. 



are the same as those of a M and in the particular 
tetrad frame and coordinate system we are using, but 
would be different in a different tetrad frame or a different 
coordinate system. 

In general, the vierbein em* 1 and its inverse e m M 
provide the means of transforming the components p M 
or of any arbitrary 4-vector between the coordinate 



frame and the tretrad frame 



Pn 



p 



(39) 



Indices on vectors p m and p m in the tetrad frame are 
raised and lowered with the Minkowski metric rj mn , 
whereas indices on vectors p^ and p M in the coordinate 
frame are raised and lowered with the coordinate metric 



9V- 



As a particular case of equations (|39H . it is true that 



(40) 



which reduces to the asserted definitions l|38ll thanks to 
the orthogonality of and /3 M . If the coordinate system 
or tetrad frame is changed, then the vierbein change 
accordingly, and a m and (3 m change in accordance with 
equations (jIU|l . 

The components u m of the 4-velocity of a particle 
relative to the the tetrad frame are related to the 
components u M in the coordinate frame by u m = e m M u M , 



(38) or explicitly 



.0 



U = V 

u l = v i - fta^ 



1,2,3) 



(41) 
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Equations (|41|l say that if in an interval of proper time 
St the particle moves a coordinate distance Sx^ — v^St, 
then relative to the tetrad frame, that is, relative to the 
locally inertial frame of an observer who is comoving with 
the infalling river, the particle moves a proper distance 



multipliers \p,v in favor of v° = 1, 



S£_ m = e m Jx' 1 = Sx m - f3 m ajx' j - 



One recognizes the right hand side of equation Ij42(l 
as having the same form as a factor of the Doran- 
Cartesian metric <|30[) . The temporal displacement S£° 
of the particle in the tetrad frame is the Galilean time 
change Stg, as in the spherical case. However, the proper 
spatial displacement S£ % of the particle in the tetrad 
frame differs from the displacement Sx % in the coordinate 
frame not by the Galilean distance f3 z Stg that the river 
moves in time Stg, as in the spherical case, but rather 
by ^a^Sx^ — (3 l (St& — asin 2 6>(50ff). The extra part 
—@ l a sm 2 9 Scftg arises from the spatial shear in the metric, 
illustrated in Figure |2 



D. Horizons 

It is now possible to see how the position of horizons 
is set by [3 = ±1, as earlier asserted, equation ijT7|) . It 
follows from the previous paragraph that the effective 
velocity of the river, from the point of view of an object 
in the river, depends on the state of motion of the object. 
The effective river velocity is f3 z a fl dx fl /dtg, which differs 
from f3 l by the factor a^v^/v — a^dx^/dtg = 1 — 
a sin 2 9 d(pg / dtg . Irrespective of this factor, the effective 
river velocity is always pointed radially (along lines of 
constant 9 and <ps) inward along the direction of j3 l . If 
we restrict temporarily to considering only objects with 
a given fixed value of a^v^/v , then such objects can 
escape outward only if their radial velocity 



dr 

~dx^' 



(43) 



exceeds zero. To determine the position of the horizon, 
we may thus first solve the slightly more general 
problem of maximizing the radial velocity v r subject to 
constraints on v° (which can be set to 1 without loss of 
generality), a^v^ 1 , and v^v^ 1 , the last constraint coming 
from the fact that the 4-velocity must be time-like or 
light-like, requiring v^v^ < 0. Equivalently, we can 
minimize v^v^ subject to constraints on v , a^v^, and 
v r . This implies that the 4-velocity must satisfy 



= A<5° + pa^ + v— 



(44) 



where A, /i, v are Lagrange multipliers, whose values are 
determined by fixing any three of the four quantities 

Not surprisingly, the largest 



v", v 



and v^v^. 



value of v r at fixed v° and a^v^ occurs when the 4- 
velocity is light-like, v^v^ — 0. Eliminating the Lagrange 



= 0, yields 



OL..V 



[R±asm9(l- (3 2 ) 1 ' 2 ] 
R{p 2 + /3 2 a 2 sin 2 (9) 



0, and 



(45) 



(42) 

which has a real solution provided that /3 2 < 1, with 
a^v^/v — p 2 /R 2 at ft 2 = 1. The position of the horizon 
is thus set by 1 = 1, as claimed: if (3 2 < 1, then there 
are geodesies on which a particle can escape, v r > 0; if 
on the other hand 1 > 1, then all geodesies are trapped, 
and an object is compelled to fall inward (or outward, in 
the case of a white hole). 

The 4-velocity of a photon that just holds steady on 
the horizon, a member of the outgoing principal null 
congruence, satisfies = (p 2 /R 2 ) dr/dx^, and is 



, ay ax 
^ = [ 1 — 



(46) 



Interestingly, the contravariant components u M of this 
4-velocity coincide, modulo a minus sign, with the 
covariant components of the azimuthal vector, 
equation Relative to the river frame, the horizon 

rotates right-handedly with angular velocity 



dt a 



a 

'R 2 



(47) 



which is also the angular velocity of the horizon perceived 
by an observer at rest at infinity. 



E. Equations of motion in the tetrad formalism 

Our aim in this subsection is to derive equations 
of motion for objects relative to the inflowing river of 
space. For clarity and pedagogy, we start from basic 
principles to derive the equations of motion H61(l of 4- 
vectors in the tetrad frame. Having derived the equations 
of motion l|61|l . we will describe what these equations 
mean physically. In the next subsection, £1111 Fl we will go 
on to apply these equations to the particular case of black 
holes, where the vierbein are given by equation (|36p . 

Let p be an arbitrary 4-vector. The 4-vector p = 
V m lm = P^SV is an invariant object, independent of 
the choice of tetrad or coordinate system. According to 
the Principle of Equivalence, an unaccelerated 4-vector p 
remains at rest in its own free-fall frame, meaning that 
its derivative with respect to its own proper time r is 
zero in its own frame 



dp 

~dT~ 







(48) 



If the 4-vector p is experiencing an acceleration in its 
own frame (perhaps because of an electromagnetic field, 
or perhaps because of rockets being fired), then the zero 
on the right hand side of equation l|48|l should be replaced 
by an appropriate invariant acceleration 4-vector. Here 
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we set any such acceleration to zero, recognizing that an 
acceleration could be reinstated if desired at the end of 
the calculation. Since p is invariant, equation (j48(l must 
be true in all frames. In the tetrad frame, this implies 



Irn , T , V u • 

dr cIt 

The proper time derivative d/dr can be written 



(49) 



d d d 

-r=v u —= u n e n v — = u n d n (50) 
dr dx v dx v K ' 



where the directed derivative d n is defined by 

„ d 



On 



dx u 



(51) 



The derivative d n defined by equation (|51|l is independent 
of the choice of coordinates x v , as suggested by the 
absence of any greek index. The derivative may be 
written d n = ~f n ■ d where d = g v djdx v and g v = 
9 v>1 9v i which shows that d n is a directed derivative along 
7„, the dot product of the vector j n with the vector 
derivative d = g v djdx v , a coordinate-independent 
object. In other words, d n constitute the tetrad frame 
components of the invariant 4-vector derivative d = 
7™ d n = g v d/dx v '. Unlike the partial derivatives djdx v , 
the directed derivatives d n do not commute. In terms of 
the vierbein derivatives d^ mn defined by 



4mn = rjki e a e 



dx v 



(52) 



the commutator [dk , d m ] of two directed derivatives is 

[^fc?^m] — fkm &n > fkmn — d nm k ^nkm • (^3) 

The fkmn are the structure coefficients of the commuta- 
tors of directed derivatives. 

Introduce the tetrad frame connection coefficients 
T mn , also known as the Ricci rotation coefficients, defined 
by 



(54) 



In terms of the vierbein e m ^ and basis vectors 
the tetrad frame connection coefficients with all indices 
lowered, Tkmn = VkiF l mn , are, from equation (|54|l . 

± kmn — 7fc ' 6Vi7m — ^k Qk ' &n 7{ ~ ■ 

ox v 



The usual coordinate frame connection coefficients, the 
Christoffcl symbols Y KUV = g K \T^, are defined by 

- 1" a 



dx v 



(56) 



Equations H55|l and (|56l) imply that the tetrad frame 
connection coefficients Tk mn are related to the Christoffcl 
symbols T Knl , by 



r. 



1. K 



(57) 



The definition lf51|l and the fact that dnijk • 7m) = 
dnVkm = implies that the tetrad frame connection 
coefficients Tkmn are antisymmetric in their first two 
indices, 



Tkmn T 



mkn • 



(58) 



The tangent vectors g^ can be regarded as coordinate 
derivatives of the invariant 4-vector interval dx = g^dx^, 
that is, g M = dx/dx^, and the commutativity of partial 
derivatives, dg^jdx u = d 2 x/dx"dx tl = d 2 x/dx fI dx' y = 
dg^/dx^, implies that the Christoffel symbols are 
symmetric in their last two indices, 

1%, = (59) 

which is the usual no-torsion condition of general rela- 
tivity. Combining equation l|57|) with the antisymmetry 
relation l|58l) and the no-torsion condition l|59l) yields an 
expression for the tetrad frame connection coefficients 
entirely in terms of the vierbein derivatives dkmn 



kmn — 2 (.dkn 



(60) 

From equations l(4^|) , (JEUJ and (|S"P it follows that the 
equations of motion for the tetrad components p k of an 
unaccelerated 4-vector p = p k ~fk are 



dp k 
~dV 



r k mn u n p m = Q 



(61) 



The physical significance of the equations of motion l|61|) 
is as follows. The tetrad 7 m defines a set of locally 
inertial frames throughout spacetime. In the present 
case, these locally inertial frames have been constructed 
so that an observer who free-falls from zero velocity 
at infinity finds their own frame aligned at all times 
with the tetrad frame. But in general another observer 
who free-falls along a different geodesic will find their 
own locally inertial frame becoming misaligned with the 
tetrad frame. Equation l|61|l expresses this misalignment 
of locally inertial frames. Because the misalignment 
is between locally inertial frames, it is a Lorentz 
transformation. This Lorentz transformation is encoded 
in the connections T k 7ln . Specifically, if a 4-vector p k is 
transported in free-fall by an infinitesimal distance 5£ n = 
u ,1 St relative to the tetrad frame j n , then the 4-vector 
experiences an infinitesimal Lorentz transformation p k — > 
pk _ S^ n T k nn p m . In other words, the connection 
coefficients T k nn for each final index n is the generator 
of a Lorentz transformation. 

The antisymmetry of the tetrad frame connection coef- 
ficient with respect to its first two indices, equation 15811 . 
expresses mathematically the fact that T k nn for each 
given n is the generator of a Lorentz transformation. 
Components of T k nn in which one of the first two indices 
k or m is (time) generate Lorentz boosts. Components 
of r^ n in which both of the first two indices k and m are 
1, 2, or 3 (space) generate spatial rotations. 
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F. The flat background 



The previous subsection, Alii El considered the equa- 
tions of motion in the tetrad formalism in the general 
case. We now particularize to the case at hand, that 
of rotating black holes, where the vierbein are given by 
equation 13tifl . In this subsection we see how the Doran- 
Cartesian coordinate system emerges as the coordinate 
system of a fiat background. In the next subsection, 
illll Gl we will see how the connection coefficients are 
expressible as the flat space gradient of a river field. In 
3III II we will revisit the notion of the flat background 
and what it means. 

Explicit computation of the connection coefficients, 
equation (JSDJ, from the vierbein of equation reveals 
that the sea of terms nonlinear in the vierbeins undergo 
a remarkable cancellation (this is not just the Jacobi 
identity at work) leaving only terms linear in the 
vierbeins e m A . In other words, the connection coefficients 
reduce to the same expression as ()60|l . but with the 
s l \ e n u factors in equation i|52|) for dfe mn replaced by 
Kronecker deltas S l x 8^ 



dkr, 



n Ox" 



(62) 



The fact that the derivative e n " d/dx v in equation i|52|) 
gets replaced by S^d/dx" in equation (|fl2*)l motivates 
introducing a new set of flat space coordinates x n , with 
latin indices, with the defining property that in the 
particular coordinate and tetrad frame that we are using 



Notwithstanding the index notation, neither x n nor x v 
is a 4-vector either under local gauge transformations 
of the tetrad or under general transformations of the 
coordinates (only the differentials dx n and dx v are 4- 
vectors), so equation l|66l) cannot be interpreted as 
a covariant equation relating the coordinates x n and 
i", even if the Kronecker delta is replaced according 
to equations 1)65(1. Rather, equation l|66|l should be 
interpreted as true in the particular coordinate and 
tetrad frame that we are using. Equation i|66f) can be 
regarded as defining the flat space coordinates x n : they 
are numerically the same as the curved space coordinates 
x v of the Doran-Cartesian metric H30|) . but reincarnated 
as coordinates x n of a flat space with a Minkowski 
metric. The Dorani^ coordinate system thus emerges 
as a rather special one, providing the coordinates of the 
flat background through which the river of space flows in 
rotating black holes. 

The flat spacetime coordinates x n are not the same 
as the locally inertial coordinates £™ attached to the 
tetrad -y n at each point of spacetime. The locally inertial 
differentials g?£™ are related to the coordinate differentials 
dx v by 



dC = e n v dx v 



(67) 



which differs from corresponding relation l|64|) between 
dx n and dx v . 



G. The river field 



dx 



~ n dx" 



(63) 



The invariant relation dx n d/dx n = dx v djdx u then 
implies that the flat space differentials dx n are related 
to the coordinate differentials dx v by 



dx n = <5™ dx" 



(64) 



It should be emphasized that the relations l|63|) and 164|) 
should be interpreted as being true only in the particular 
tetrad and coordinate frame that we are using. If the 
tetrad frame is subjected to a local gauge transformation 
(i.e. a Lorentz transformation that varies from place to 
place) that rotates the locally inertial coordinates at each 
point by £ n — > £' n , and if the coordinate system is 
subjected to a general coordinate transformation x u — > 
x' v , then the Kronecker deltas in equations l|63|l and l|64|) 
should be replaced by 



(65) 



In the particular tetrad and coordinate frame that we are 
using, integrating the relation (|64|) arbitrarily through 
space yields (the constant of integration being set to zero 
without loss of generality) 

x n = 8™ x v . (66) 



The vectors a m and /3 m can be regarded as functions 
of the flat space coordinates x n , and the replacement 
of the vierbein derivatives dkmm equation l]62p. in the 
connection coefficients can be written 



'■km n 



da m Pk 
dx n 



(68) 



The connection coefficients, equation (|60f) . are then given 
by flat space derivatives of a m and /3 m 



da m (3 k da k [3 m da m f3 n da n (3 k 



dx 11 
da„(3r, 
dx k 



dx r 
dx m 



dx k 



dx n 



(69) 



The connection coefficients with zero final index n = 
are all identically zero, rfc m o = 0, and taking the 
spatial curl of Tkmn on the n index yields another sea 
of terms which again undergo a remarkable cancellation 
to nothing 



krrin 



o 



(70) 



for all i,fc,m. This demonstrates that the connection 
coefficients Yk mn must be expressible as (minus) the flat 
space gradient d/dx" of an object u>k m , which we call the 
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FIG. 3: (Color online) The velocity and twist fields for an 
uncharged (Kerr) black hole with angular momentum per 
unit mass a = 0.96. The arrowed lines show the magnitude 
and direction of the river velocity, while the unarrowed lines 
emerging from the arrowed lines show the magnitude and axis 
of the river twist. The confocal ellipses show the outer and 
inner horizons, and the large dots at the foci of the ellipses 
indicate the ring singularity. In the vacuum Kerr solution, the 
river velocity goes to zero at the horizontal disc bounded by 
the ring singularity, then turns around and rebounds through 
a white hole into a new universe. 



while the angle and axis of the river twist are given by 
the rotation vector 

= \ e lkm uj km = e ikm a k f3 m + C (*, k, m = 1, 2, 3) . 

(75) 

Like the velocity vector /3 i; the twist vector /j? at each 
point lies in the plane of constant free-fall azimuthal angle 
(j>g, since it is a sum of two vectors e lkm a k p m and Q both 
of which are orthogonal to the azimuthal vector a k ■ 

Figure illustrates the velocity and twist fields Pi 
and fii for an uncharged black hole with specific angular 
momentum a — 0.96. 

Another familiar bivector is the electromagnetic field 
tensor F km , and it can be useful to think of the river 
field bivector io km in those terms. The velocity vector 
Pi is the analog of the electric field vector Ei, while the 
twist vector fa is the analog of the magnetic field vector 
Bi. The analogy extends to the fact that, like a static 
electric field, the velocity vector pi is the gradient of a 
potential ip, 

A = -|t, ^^-j^Pdr. (76) 

However, unlike a magnetic field, the twist vector \i l is 
not pure curl, although curiously (J, 1 + Q l is pure curl, 
having zero divergence, d(fj, z + (*)/dx % = 0. 



river field since it encapsulates all the properties of the 
river in the river model: 



duj kr , 

dx n 



(71) 



The river field 0Jk m is a bivector^L^i, inheriting from 
rfcmn the property of being antisymmetric in km. That 
the connection coefficient T kmn is the flat space gradient 
of the river field lies at the heart of the river model as a 
description of black holes. After some manipulation we 
find the desired bivector river field to be 



^>km — CtkPm — &mPk + EQkmi C 



where the vector £ l is 



C = (o,o,o,c) , C = 



Pdr 



(72) 



(73) 



which points vertically upward along the rotation axis of 
the black hole. 

The river field io km given by equation (|72|l inherits from 
the connection coefficient T kmn its Lorentz structure. 
The river field defines a velocity and a rotation, or twist, 
at each point of the black hole geometry. Components of 
bJkm i n which one of the indices k or m is (time) define 
a velocity, while components in which both indices k and 
m are 1,2,3 (space) define a spatial rotation, or twist. 
The velocity is just the river velocity p m 



H. Motion of objects in the river 

We are now ready to demonstrate a fundamental 
feature of the river model for stationary rotating black 
holes, that as an object moves through the river of space, 
it is Lorentz boosted and rotated by the tidal gradients 
in the velocity and twist fields of the river. 

It follows from inserting the connection coefficients 
Tkmn from equation 171|) into the equation of motion 161|) 
that the equation of motion of an unaccelerated 4-vector 



• k in the river frame is 



~dV 



dui k 7 



-u n p m 



(77) 



The equation of motion l|77|l can be interpreted as follows. 
In an infinitesimal interval 5t of proper time, a particle 
moves a distance S£ n = u n Sr relative to the infalling river 
of space. As a result of its motion through the river, the 
particle experiences a tidal change 



duj k 



-6C 



(78) 



in the river field, which generalizes equation 119|) for 
spherical black holes. The tidal change Sui k m in the river 
field is an infinitesimal Lorentz transformation, and it 
induces a Lorentz boost and rotation in the 4-vector p k 



<^>0r, 



Pn 



(74) 



P 



p k + 5uj k m p r ' 



(79) 
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FIG. 4: (Color online) Two sample sets of geodesies in 
the equatorial x-y plane (0 — 7r/2) of an uncharged (Kerr) 
black hole with angular momentum per unit mass a = 0.96, 
plotted in Doran-Cartesian coordinates, frill Bl Each set of 
geodesies shows a central point (large dot) and associated 
locally inertial axes (crossed thick lines), surrounded by a 
set of points that are initially uniformly spaced around, and 
initially at rest relative to, the central point, in the locally 
inertial frame of the central point. The central point and its 
attendants follow geodesies (thin lines) into the black hole, 
the ensemble becoming tidally distorted as it falls in. In the 
upper set of geodesies, the central point is comoving with the 
infalling river of space, while in the lower set of geodesies, the 
central point is initially moving radially outward, but soon 
turns around and falls in. The lower ensemble illustrates how 
the locally inertial axes attached to the central point twist 
as the ensemble falls in, the twist acting so as to keep the 
frame comoving with the geodesic motion of points in a small 
neighborhood of the central point. 



Equations (|78|l and (|79|l reproduce the equations of 
motion ijTTj) . 

Figure 01 shows two ensembles of geodesies computed 
using the equation of motion (|77|l . Each ensemble 
consists of a central point and associated tetrad axes 
surrounded by a set of points that are initially uniformly 
spaced around, and initially at rest relative to, the central 
point, in the locally inertial frame of the central point. In 
the upper ensemble, the central point is comoving with 
the infalling river of space, while in the lower ensemble 
the central point is initially moving radially outward, but 
soon turns around and falls inward. The tetrad axes are 
skewed because the spatial metric is sheared (compare 
Figure In the lower ensemble the tetrad axes are 
Lorentz-contracted in the radial direction because of the 
initial outward motion of the ensemble relative to the 
infalling river. The ensembles of points become tidally 
distorted as they fall into the black hole. If the locally 
inertial coordinates of a tetrad axis are denoted S^ k , 



then the tetrad axis evolves according to the equation 
of motion (|77J) with p k — > 5£ k , 



d5£ k du k 



dr 



dx r - 



l -u n SC 



(80) 



Similarly, the tetrad 4-velocity u k of each point in the 
ensemble evolves according to the equation of motion 177|) 
with p k — > u k , 



du k 
~dr~ 



duo k 



dx n 



m u n u m 



(81) 



Each point surrounding the central point is initially at 
rest relative to the central point, in the latter's locally 
inertial frame. This requires that the covariant difference 
in tetrad 4-velocities between each point and the central 
point initially vanishes, which requires that the difference 
Su in the tetrad 4-velocity of a point separated from the 
central point by locally inertial separation 5£ k initially 
satisfies, to linear order in the separation 5£ , 



5u k 



dx r - 



(82) 



The difference Su k in equation (|82|l is to be understood 
as the tetrad 4-velocity of a point evaluated in the tetrad 
frame at that point, minus the tetrad 4-velocity of the 
central point evaluated in the tetrad frame at the central 
point. Notice that the indices on S£ n and u m on the 
right hand side of equation (|52*|l are swapped compared 
to those on the right hand side of equation l|80|) . In 
equation (|80|l the axis <5£ fe is transported along the 4- 
velocity it", whereas in equation l)82|) the 4-velocity u k is 
transported along the axis 6£ n . 

The lower ensemble of points in Figure 0] illustrates 
the twist in the locally inertial frame that develops as 
the ensemble moves through the river of space. The twist 
acts so as to keep the locally inertial frame comoving with 
the geodesic motion of points in a small neighborhood of 
the frame. 

Equation 182fl is true initially, when the ensemble 
of points are at rest relative to each other, satisfying 
D5^ k / Dt = 0. The more general form of equation 182|l . 
valid when the points are in relative motion is, to linear 
order in the separation 5£ k (here we revert, eq. (I71II . to 
the more familiar notation for the connections r^ n ), 



5u k + r k m jCu m = 



DS£_ k 
Dt 



(83) 



or equivalcntly 
Su k = 



k _ dse 

dr 



lyK _ yk \ Sjtn m 
,A ram L mn) °S u 



(84) 



Variation of the equation of motion l|81|l for u k gives 



dSu k 

~dT 



+ diT k m Ji l u n u m + T k nn (u n Su m + u m Su n ) 



. 

(85) 
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Inserting the expression (|84|l for 5u k into this equa- 
tion 185|) yields the familiar equation of geodesic 
deviation 

= R lmn k u m u n 5£! (86) 
where Rkimn is the Riemann curvature tensor 

Rklmn — &l^nmk ^fcfnml ~t~ L m fJ- jnk ^ mk^ 3 n ^ 

+ ( r /fe _ r L) r «mj ' ( 87 ) 

I. The flat background revisited 

Now that we have completed the formalism of the 
river model, it is useful to revisit the question of the 
flat background, flIII Fl through which the river of space 
flows and twists into a rotating black hole. What exactly 
does flatness mean in this context? 

The crucial equation is equation Ij71|l . which states 
that the connection coefficient is given by the flat 
space gradient of the river field. The fact that the 
gradient is an ordinary partial derivative with respect 
to Doran-Cartesian coordinates is what makes the 
background flat, and in a sense that is all there is to 
it. Equation l|71[) acquires physical significance because 
it propagates through to the equation of motion (|77ll of 
objects swimming in the river. The equation of motion 
paints the physical picture of objects moving in the river 
being Lorentz boosted and rotated by the flat space tidal 
gradients in the velocity and twist components of the 
river field. 

The statement that the background spacetime in the 
river model is flat is not a statement about the metric 
g^ v being flat. Rulers and clocks swimming in the river of 
space measure not distances and times in the background 
space, but rather distances and times relative to the 
tidally twisting and stretching river. The presence of 
tides is the signature of curvature, so it makes sense that 
the metric measured by rulers and clocks is not flat. 

It is to be emphasized that the flat background has no 
physically observable meaning. It is simply a fictitious 
construct that somehow emerges from the mathematics. 

IV. SUMMARY 

In this paper we have presented a way to conceptualize 
stationary black holes, which we call the river model. 
The river model offers a mental picture of black holes 
which is intuitively appealing, and whose basic elements 
are simple enough that they can be grasped by non- 
experts. In the river model, space itself flows like a river 
through a flat background, while objects move through 
the river according to the rules of special relativity. For 



a Schwarzschild (non-rotating, uncharged) black hole, 
the river falls radially inward at the Newtonian escape 
velocity, hitting the speed of light at the horizon. Inside 
horizons, the river of space moves faster than light, 
carrying everything with it. 

We have presented the details that place the river 
model on a sound mathematical basis. We have shown 
that the river model works for any stationary black 
hole, rotating as well as non-rotating, charged as well 
as uncharged. The Doran^S, coordinate system provides 
the coordinates of the flat background through which the 
river of space flows into the black hole. 

The extension of the river model to rotating black holes 
proves to be both surprising and pretty. Contrary to 
expectation, the river does not spiral into a rotating black 
hole: the azimuthal component of the river velocity is 
zero. Instead, the river has at each point not only a 
velocity, but also a rotation, or twist. The river is thus 
a Lorentz river, characterized by all six generators of the 
Lorentz group. As an object moves through the river of 
space, it is Lorentz boosted by changes in the velocity 
of the river along its path, and rotated by changes in 
the twist of the river. Equation (|72l) gives an explicit 
expression for the river field, a six-component bivector 
field that specifies the velocity and twist of the river at 
each point of the black hole geometry. 

The tidal boosts and twists experienced by an object 
in the river induce a curvature in the spacetime measured 
by the object, causing the metric to be non-flat. Changes 
in the river velocity rotate between space and time axes, 
while changes in the river twist rotate between two 
spatial axes. For a spherical black hole, the river has 
zero twist, so objects experience no spatial rotation, with 
the consequence that the metric, the Gullstrand-Painleve 
metric, is flat along spatial hypersurfaces at constant 
time, dtff = 0. For a rotating black hole, the river has 
a finite twist, and the metric is not flat along spatial 
hypersurfaces. 
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APPENDIX A: PROJECT: THE RIVER MODEL 
OF BLACK HOLES 

The project below has been field-tested and refined 
over a period of several years in undergraduate classes 
on relativity and black holes at both lower-division non- 
science-major and upper-division science-major levels. It 
was designed as a 45-minute "in class group project" , in 
which students would split into groups of 3 or 4, and 
by arguing with each other would arrive at consensus 
answers to a series of concept questions. At the end of the 
project each group would submit its answers for grade. 

CONCEPT QUESTIONS 

According to the river model of black holes, the 
behavior of objects near black holes is precisely as if 
space were falling like a river into the black hole. For 
spherical black holes, this model was discovered in 1921 
by the German Nobel prizewinner Allvar Gullstrand and 
independently by the French mathematician Painleve. In 
the model, space falls inward at the Newtonian escape 
velocity v = y/2GM/r. The infall velocity is less than 
the speed of light c outside the horizon, equals the speed 
of light c at the horizon, and exceeds the speed of light c 
inside the horizon. 

What does the river model predict for the answers to 
the questions below? [For freshman non-science majors, 
use only the unstarred questions. For more advanced, 
science-major students, use all questions, and drop or 
abbreviate the hints.] 

*1. What radius does the river model predict for the 
horizon of a black hole? 

2. Suppose that you are a light beam (therefore 
moving at the speed of light) exactly at the horizon. 
What would happen to you if were pointed directly 
outward? [Do you fall in? Do you move out? Do 
you move sideways?] What would happen to you if 
you were pointed mostly but not exactly outward? 

3. In what way, if any, does this behavior differ from 
the predictions of the corpuscular theory of light, 



which in the hands of John Michell in 1784 gave 
the "correct" result for the radius of the horizon? 
[In the corpuscular theory of light, a corpuscle 
of light is emitted at the speed of light, and 
thereafter behaves much like a massive particle: it 
flies outward, and it either goes to infinity or turns 
around and comes back depending on whether its 
initial velocity, the speed of light, is more or less 
than the escape velocity] 

4. Suppose that you are a light beam orbiting the 
black hole in a circular orbit. On this orbit, the 
so-called "photon sphere" , are you at the horizon, 
inside the horizon, or outside the horizon? Justify 
your answer. 

5. Make a connection between the appearance of the 
sky if you hover just above the horizon of a black 
hole, and special relativistic beaming. [How does a 
scene appear if you move through it at very close 
to the speed of light?] 

6. Qualitatively, what would the river model predict 
for the tidal forces experienced by an infalling 
observer? [First, the tidal force in the vertical 
direction. Think about the fact that the river is 
accelerating inward. Next, the tidal force in the 
horizontal direction. Think about the fact that the 
river is converging (getting narrower) as it flows 
inward.] 

*7. How does the river model account for redshifting 
and freezing at the horizon? 

*8. Given that one of the fundamental propositions of 
Special and General Relativity is that spacetime 
has no absolute existence, what does it mean to 
say that space is falling into a black hole? 

*9. In the river model, the flow of space accelerates 
inward to the black hole. If the river were moving 
uniformly instead of accelerating, would there be 
any gravity? 

ANSWERS 

1. The river velocity equals the speed of light 

1/2 

when (2GM/r s y^ = c, which rearranges to an 
expression for the radius of the event horizon, the 
Schwarzschild radius r s , 

2GM 
r s = — — . 

<r 

2. If you were a light beam pointed directly outward 
at the horizon, then you would hang forever at the 
horizon, your outward motion at the speed of light 
being exactly canceled by the inward motion of the 
river of space at the speed of light. If you were 
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a light beam not exactly pointed outward, then 
the outward component of your velocity would be 
a bit less than the speed of light, since part of your 
velocity would be sideways. The inflow of space 
would then carry you into the black hole. 

3. Whereas in general relativity an outwardly pointed 
light beam at the horizon hangs there motionless 
for ever, in the classical corpuscular theory the light 
never remains at rest. The light either keeps going 
outward for ever (if its velocity exceeds the escape 
velocity), or else it turns around and comes back. 
It is true that the light is motionless at the instant 
of turnaround, but otherwise the light is always 
moving. 

Another difference is that in general relativity the 
question of whether a light beam can escape from 
a point just above the horizon depends on the 
direction in which the light beam is pointed. If 
the light beam is pointed directly outward, then it 
will escape, but if it is pointed somewhat sideways, 
then it will fall into the black hole. In the 
classical corpuscular theory, by contrast, whether a 
corpuscle escapes from a given point depends only 
on whether its velocity exceeds the escape velocity, 
not on the direction in which it is pointed. 

4. You cannot be at the horizon, because if you had 
any sideways motion, which you must because you 
are in circular orbit, then the inflow of space would 
drag you into the black hole. And you cannot 
be inside the horizon, because the inflow of space 
would again drag you inwards. Therefore you must 
be in circular orbit somewhere outside the horizon. 
For a Schwarzschild black hole, the radius of the 
photon sphere turns out to be 1.5 Schwarzschild 
radii. 

5. If you move through a scene at very close to 
the speed of light, then the scene ahead of 
you, in the direction you arc moving, appears 
concentrated, brightened, and blueshifted. If you 
hover just above the horizon of a black hole, 
then according to the river model you must be 
moving very rapidly through the inflowing river 
of space. Consequently the view above you must 
appear concentrated, brightened, and blueshifted. 



It should be emphasized that hovering just above 
the horizon of a black hole is an unnatural and 
wasteful thing to do. In reality, you would surely 
"go with the flow" of space. If you free-fall into 
a black hole, then you do not see the sky highly 
concentrated above you. 

6. Since the river is accelerating inwards, the velocity 
of the river is faster at your feet than at your head 
(presuming that you are upright, so that your feet 
are closer to the black hole than your head). The 
difference in river velocity means that you feel a 
tidal force in the vertical direction, pulling your 
feet away from your head. 

In the horizontal direction, the river is converging 
spherically towards the black hole, so you feel 
tidally squashed in the horizontal direction. 

7. Just above the horizon, a photon battling against 
the inrushing torrent of space takes a long time to 
get to an outside observer. As the emitter gets 
closer to the horizon, it takes longer and longer 
for the photon to get out, until at the horizon it 
takes an infinite time for a photon to lift off the 
horizon. Thus as an object approaches the horizon, 
it appears to an outside observer slower and slower, 
thus more and more redshiftcd. Asymptotically, 
the object appears to freeze on the horizon, and 
the redshift tends to infinity. 

8. The river model consists of a set of coordinates 
(the background) and a set of locally inertial frames 
that flow through those coordinates (the river that 
flows through the background). Attaching a set of 
coordinates and a set of locally inertial frames does 
not make the spacetime absolute. 

9. According to the Principle of Equivalence, a 
gravitating frame is equivalent to an accelerating 
frame, so if there is no acceleration, then there is no 
gravity. However, if the river is falling at constant 
velocity in the vertical direction but still converging 
horizontally because of the spherical convergence of 
the flow, then you will feel a tidal squashing in the 
horizontal direction, so there must be a gravity. 



